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Z ' Abstract: 

^ ! The effect of bulk viscisity on the evolution of the homogeneous 

and isotropic cosmological models is considered. Solutions are 
found, with a barotropic equation of state, and a viscosity coef- 
ficient that is proportional to a power of the energy density of 

t^^ . the universe. For flat space, power law expansions, related to 

extended inflation are found as well as exponential solutions, re- 
lated to old inflation; also a solution with expansion that is an 
exponential of an exponential of the time is found. 
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1 Introduction 

The effect of negative pressure and bulk viscosity on the evolution of homo- 
geneous isotropic cosmological models have been considered recently by Wolf 
[11]. He considered a viscosity coefficient that is proportional to the energy 
density of the matter in the universe and fiat topology of the three space. In 
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this paper we extend his consideration to the case when the viscosity coeffi- 
cient is proportional to a power of the energy density and arbitrary curvature 
of the three space. 

We want to consider Eintein's field equations for the isotropic and homo- 
geneous line element, 



ds^ = dr - a\t) 



1 — fcr^ 



+ r'{de' + sin' 



(1) 



and with the material substratum corresponding to a fluid with bulk viscosity, 
whose energy momentum tensor is given by 



T^u = {p + P- COWu, + p- c%., e = u%. (2) 



They are 



d 2 k 9)71 pG 
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(3) 



- + (-)^ + 4 = -8^G(P-3C-), (4) 

a a a^ a 

eliminating the energy density from the above equations we obtain, 

- + (-r + 4 + (-)' = -8-G(^-3C-). (5) 

a a a"^ a a 

In order to continue we need some equation of state in the form p = p{p) and 
C = C(p) that we assume of the following form, 

P = ep, C = Cop". (6) 

The last equation of state has been used before. If n=l we have a relativis- 
tic fluid, n=3/2 corresponds to a string dominated universe in the sense of 



Turok [0] . Some other values of n arise when quantum effects are considered. 
With the above equations of state Eq.(^) reduces to, 



2ad + A{d^ + k)+B^{a' + k)'' = 0, (7) 

(Jj 



where 



A = 3e + 1, B = -(87r)(i-")3(i+")Co- (8) 

After solving the above equation we use Eq.(^) to calculate the corre- 
sponding energy density of the models. 

To solve the differential equation it is useful to define the following de- 
pendent variable, 

y = y{a) = a" + k (9) 

and now our differential equation is 

ay' + Ay + Ea^^"'") {y - k) ^/'y" = 0. (10) 

In the next sections we solve exactly this equation for some values of A, B 
and n. 



2 Flat Space Solutions 

For the case of flat space (k=0), Eq. (p!oD reduces to 

2aa + Aa^ + 5a^~2n^i+2n _ q^ ^^^^ 



ay' + Ay + Ba^^-^"''>y^/^+'' = 0. (12) 

2.1 n=l/2 

For n = 1/2 the above equation is linear, 

ay' + Cy = 0, C = A + B = 3e + 1 - {8nf/^h'^^^^\o, (13) 

and we can find a(t) explicitely. For C = — 2 we have, 

a = aoe^*, for C = -2, (Stt) (1/2)3(3/2)^^ ^ 3^^ ^ ^^^ ^^4^, 

and the corresponding energy density is 

and for C 7^ —2 the solution is 

a = [ki + k2t]^ , ioT C j^ -2 (16) 

with the density 

where /Cj, Oq, and H are constants of integrations. The first of the above 
solutions is an old inflationary solution, the second one for any (^0 7^ will 
correspond to the new infationary solutions. 



2.2 717^1/2,^7^-2,(6 7^-1) 



For n 7^ 1/2 Eq.(|ToD is of the Bernoulli type and the solution for A ^ —2, (e ^ 
-1), is 

y = [yia(^"-i)^/^ - lf2^'''"]'^' ^ ^ -2 (^ ^ "l)" (l^) 

For n=l it is possible to do the second integration for arbitrary value of 
the integration constant yi, the result is, 

2yia'-^ + 9^0 ln(a) = ±3(e + 1) (t + to) (19) 

with the density given by 

3 . 3(.+i) 9?7o T_2 /„^s 

This solution has the exponential behaviour for earlier times, with Hub- 
ble parameter H = {e + l)/3^0; and a power law expansion at late times 
with exponent equal to 2/3 (e + 1). This late behaviour corresponds to the 
standard cosmology without viscosity. Therefore this family of solutions pro- 
vides a natural end to the inflationary era with exponential expansion with 
a transition to a standard cosmological era if e > —1/3 or to an extended 
inflatinary period if e < — 1/3. 

To proceed for n 7^ 1 we have to make the choice yi = 0, and we have 

y = a' = [-^a^-^"]!^, A j^ -2 (e ^ -1) (21) 

from where the expansion factor can be obtained explicitely, 

a(t)=aoe^*, ^ = t^]^ (22) 



with the corresponding energy density 






2.3 n^^ 1/2,^= -2, (e= -1) 

In case our fluid is of the vacuum type, Eq. ([l0|) can be integrated, 

2/ = a' = a2[yi + ^ _^lna]T^, A = -2 (e = -1). (24) 

In this case we can do the second integration without any assumption 
about the integration constant yi 
For n y^ the solution is 

.2{ai+nBty 



a{t)=aoeM ^^^ _ ^^^ ], (25) 

oil , . 2n—2 , ^ 

P=^^i + ^B^)~- (26) 

For n=0 the solution is 

a{t) = aoexp[aiexp[—Bt/2]], (27) 

here Oj are the redefined constants of integration. From the above solutions 
in the case of flat space we see that for any equations of state ( arbitrary 
e and n ) we have inflation, i.e., accelerated expansion. In particular for 
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e = — 1, that gives the old inflation in the usual case without viscosity, we 
see that viscosity gives a " superinflation" , as long as ai > 0. 



3 Non-flat Solutions 

In the case of non-flat space Eq.(|TUp can be solved for n=0, and A=0 (e = 
— 1/3) the solution is 

a{t) = ao + aie~-^*/^ (29) 

We notice that the expansion factor does not depend on the value of the 
curvature of the three dimensional space, the reason for that can be seen from 
equation (7) since for n=A=0 the dependence on the curvature disappears. 
Nevertheless, the curvature is present in the energy density. We recall here 
that the equation of state corresponds to a gas of textures and also to some 
cases of cosmic strings. 
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